Capacities of lossy bosonic memory channels 
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We introduce a general model for a lossy bosonic memory channel and calculate the classical and the quantum 
capacity, proving that coherent state encoding is optimal. The use of a proper set of collective field variables al- 
lows to unravel the memory, showing that the rc-fold concatenation of the memory channel is unitarily equivalent 
to the direct product of n single-mode lossy bosonic channels. 
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PACS numbers: 03.67.Hk, 05.40.Ca, 42.50.-p, 89.70.-a 

One of the most important problem of quantum informa- 
tion theory is finding the maxima rates (i.e. capacities) at 
which quantum or classical information can be transmitted 
with vanishing error in the limit of large number of transmit- 
ted signals (channel uses) 0]]. Earlier works on the subject 
focused on models where the noise affecting the communica- 
tion is assumed to act independently and identically for each 
channel use (memoryless quantum channels). Recently, how- 
ever, an increasing attention has been devoted to correlated 
noise models (memory quantum channels), see e.g. yj] and 
Ref.s therein. Memory effects in the communication may 
arise when each transmitted signal statistically depends on 
both the corresponding and previous inputs. Such scenario 
applies when the dynamics of the communication line is char- 
acterized by temporal correlations which extend on timescales 
which are longer than the times between consecutive chan- 
nel uses — a regime which can be always reached by in- 
creasing the number of transferred data per second. For in- 
stance optical fibers may show relaxation times or birefrin- 
gence fluctuations times longer than the separation between 
successive light pulses J^. Similar effects occur in solid state 
implementations of quantum hardware, where memory effects 
due to low-frequency impurity noise produce substantial de- 
phasing JfU] . Furthermore, moving from the model introduced 
in memory noise effects have also been studied in the con- 
test of many-body quantum systems by relating their proper- 
ties to the correlations of the channel environmental state |6[] 
or by studying the information flow in spin networks J3l • 

It is generally believed that memory effects should improve 
the information transfer of a communication line. However 
finding optimal encodings is rather complex and up to date 
only a limited examples have been explicitly solved 
In this paper we focus on a continuous variable model of quan- 
tum memory channels in which each channel use is described 
as an independent bosonic mode. The proposed scheme is 
characterized by two parameters which enable us to describes 
different communication scenarios ranging from memoryless 
to intersymbol interference memory 1 1 Oil - up to perfect mem- 
ory configuration ll 111 . It effectively mimics the transmission 
of quantum signals along attenuating optical fibers character- 
ized by finite relaxation times, providing the first comprehen- 
sive quantum information characterization of memory effects 
in these setups. For such model we exactly calculate the clas- 




/// 
i 



FIG. 1: Left: a single use of the memory channel (see text for de- 
tails). Right: the n-fold concatenation of the memory channel: pho- 
tons entering in the fc-th input mode ak can only emerge in the output 
ports by with k' ^ k (the channel is thus non-anticipatory). 



sical and the quantum capacity 11211 and prove that coherent 
state encoding is optimal. This is accomplished by unravel- 
ing the memory effects through a proper choice of encoding 
and decoding procedures which transform the quantum chan- 
nel into a product of independent (but not identical) quantum 
maps. If the channel environment is in the vacuum, the capaci- 
ties can then be computed by using known results on memory- 
less lossy bosonic channels jl3lll4ll which in the limit of large 
channel uses provide converging lower and upper bounds. 

Channel model:- We consider quantum channels described 
by assigning a mapping of the form 



TTE[U n (p n ®<T E )Ut] 



(1) 



where p n and & n (p n ) represent, respectively, the input and 
output states of the first n channel uses, and <je is the initial 
state of channel environment E. The latter is composed by 
a memory kernel M which interacts with all inputs, and by 
a collection Ei, E2, ■ ■ ■ , E n of local environments associated 
with each individual channel use. Such interactions are de- 
scribed by the unitary U n which can be taken as a product of 
identical terms, i.e. U n = U n U n _i ■ ■ - U\ with being the 
interaction between the A;-th channel input, Ek and M . Within 
this context the channel uses will be described by an ordered 
sequence of independent bosonic modes associated with the 
input mode operators {ai, 02, • • • , a„}. Through the cou- 
pling U n they undergo a damping process that couples them 
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with the local environments E\, E2, ■ ■ ■ , E n and the memory 
kernel M (also described by a collection of mode operators 
{ex, &ii ' ' ' j e «} and mi). Memory effects arise when the 
photons lost by the fc-th channels mix with the environmen- 
tal mode ek+i of the subsequent channel use. Specifically the 
evolution of fc-th input is obtained by a concatenation of two 
beam-splitter transformations, the first with transmissivity e 
and the second with transmissivity 77 (see Fig. Q] left). In the 
Heisenberg-picture this is defined by the identities 



= y/erjm k + yjl - r/a k + y/r)(l - e) e k , 



-\Je(l- rj) m k + ^rja k - \J (I - e)(l - 77) e k , 
-Vl -em k + v^efc, 



(2) 



where m' k := ulm k U k , b k := U~la k U k , and e' k := U\e k U k 
describe the outgoing modes of the model (in particular the 
b k s are associated with the receiver signals). The resulting in- 
put/output mapping is finally obtained by a n-fold concatena- 
tion of Eq.s (|2} where, for each fc, we identify the mode m k +i 
with m' k (see Fig.Q] right). This yields a non-anticipatory IU5I1 
channel where a given input can only influence subsequent 
channel outputs (i.e. for each fc, b k depends only upon the 
dfc/'s with fc' ^ fc). The transmissivity e plays the role of 
a memory parameter. In particular the model reduces to a 
memoryless channel 111 311 for e = (the input a k only in- 
fluences the output b k ), and to a channel with perfect memory 
I n|] for e = 1 (all a k interacts only with the memory mode 
mi). Intermediate configurations are associated with values 
e s]0, 1[ and correspond to intersymbol interference channels 
where the previous input states affect the action of the chan- 
nel on the current input ifioll . Of particular interest is also the 
case 77 = where <!>„ describes a quantum shift channel ifioll . 
where each input state is replaced by the previous one. 

When dealing with memory channels, four different cases 
can be distinguished depending whom the memory mode is 
assigned to |2]. Specifically the initial and final state of the 
memory can be under the control of the sender of the mes- 
sage (A), the receiver (B) or the environment (E). The four 
possible XY setups are denoted: XY = AB (initial mem- 
ory to A and final memory to B), XY = AE, XY = EB, 
XY = EE. These different scenarios typically lead to dif- 
ferent values of the channel capacity but, at least for finite di- 
mensional system, they coincide if the channel is forgetful 10] . 
To make the notation homogeneous we thus define: ao := mi 
and b 



n+l 



if XY = AB; clq := mi and e' n+1 := m' r 



if XY = AE; e := m x and b n+1 := m' n if XY = EB; 



cq := mi and e n+1 := m' n if XY = EE. 

With the above choices the output modes of the receiver can 
then be expressed in the following compact form 



Ul a k U T 



XY 



(3) 



cumbersome]. The A k Y commute with the E^ together 
with their hermitian conjugates. Furthermore they satisfy the 
following commutation relations: 



[E k Y , E k , Y ^} 



M k x k Y, 



with 6 kk / being the Kronecker delta and M XY being a sym- 
metric, positive real matrix which satisfies the condition 1 ^ 
M XY . For example the n x n matrix M EE has elements 



M kU = S kk' - (1 - Vmin{k,k'})V^n 



k-k'\ 



with i]k := 77 + (1 - (e?7) ) 



k-l^ejl-rj) 2 



l — er} 



Analogous expres- 



sions hold for XY = AB, AE and EB which only differ 
by terms which in the limit of n — ► 00 can be neglected. In- 
deed, by varying n, the M XY form a sequence of matrices of 
increasing dimensions which (independently from XY) are 
asymptotically equivalent IU6I1 to the Toeplitz matrix M^°°^ 
of elements 



(00) 

kk' 



Skk' - (1 - r? (oo) )^ 



-\k-k'\ 



(4) 



with 77 



(00) 



liim 



• Vk = V + 



Similarly the 



asymptotic distribution of the eigenvalues t xy of M XY can 
be computed by performing the Fourier transform of the ma- 
trix Defining z := 27rfc/n and taking ?? — > 00 this 
gives the nondecreasing function 



r(z) 



e + 77 — 2^/efjcos (z/2) 
1 + erj — 2y/e?7COS (z/2) 



iz/2 



1 - ^/Wie 



iz/2 



(5) 

which is plotted in Fig. |2{c). According to the Szego theo- 
rem IU6I1 the asymptotic average of any smooth function F of 
the eigenvalues of M Y can then be computed by the formula 



^F(r(z)l 



(6) 



which is explicitly non dependent upon XY. 

Unraveling the memory:- We show that the memory effects 
can be unraveled by introducing a proper set of collective co- 
ordinates. To do so we introduce the (real) orthogonal ma- 
trix O xy which diagonalizes the matrix M XY (it exists since 
the latter is real symmetric), i.e. J2 r ,r' O k Y M X J O xy , = 
S kk 'T XY (here the r xy £ [0, 1] are intended to be arranged 
in nondecreasing order). 

Let us define the following sets of operators := 

T, k .O x k Y b k „ a fe := Z k ,0 XY A XY /^ 

£fc< o 



r XY 



'kk' 



By construction they satisfy 



canonical commutation relations, moreover it is easy to show 
that they obey the following transformations 



with A XY and E x being, respectively, field operators 
formed by linear combination of the field modes ay and e k * 
with fc' ^ fc [The explicit expressions can be easily derived 
from Eq. (O but are not reported here because they are rather 



m*kU n = Jr XY 



1 - ^ Y e k . 



(7) 



We denote by Wa, Vb, Te the canonical unitaries 017I1 that 
implement the transformations a k — > & k = W A a k WA, b k — > 
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bfe = VgbkVB and — > = T E etTE- We have shown that 
the channel <&„ is unitarily equivalent to the map 



<(pj = Tr B [c/;( Pn ® - , B )([/; 



(8) 



with er^ := T^cteTe, and where the unitary transformation 
C/'j := VaUtiI Wa ® 7b) induces the beam-splitter transfor- 
mations in (IT) 111 811 . Formally, the unitary equivalence reads 
*«(Pn) = ^n{W\p n WA)V\, i.e. we can treat the out- 
put states of <&„ as output of <&'„ by first counter-rotating the 
input p n by Wa (coding transformation) and then by rotat- 
ing the output by Va (decoding) 1 5]. Assuming then cte to be 
the vacuum state, we have cr' E = <te and the map ([8]l can be 
written as a direct product of a collection of independent lossy 
bosonic channels, i.e. 



(9) 



with being a single-mode lossy bosonic channel with ef- 
fective transmissivity t^ y . 

Classical capacity:- Equation (O suggests that we can 
compute the classical capacity of <!>„ by applying the results 
of Ref. 1 1311 on memoryless multi-mode lossy channel. To do 
so however, we have first to deal with the fact that the single- 
mode channels forming & n are not necessarily identical (in- 
deed, for finite n their transmissivities t^ y can be rather dif- 
ferent from each other). Therefore the map (O is not memo- 
ryless in the strict sense. To cope with this problem we will 
construct two collections of memoryless multi-mode channels 
which upper and lower bound the capacity of & n (and thus of 
3>„), and use the asymptotic properties of the distribution (0 
to show that for large n they converge toward the same quan- 
tity. 

First as usually done when dealing with bosonic chan- 
nels II 1 911 . we introduce a constraint on the average pho- 
ton number per mode of the inputs signals. This yields 
the inequality — J2k T r [aI a fcP n ] < N, which is pre- 
served by the encoding transformation p n — ► W A p n WA of 
Eq. dHJ due to the fact that Wa is a canonical unitary, i.e. 
i J^k ^ Y [ a k a ^APn W A\ < N. For any n, we then group 
the single-mode channels of Eq. (O in J blocks, each of size 
£ = nj J. At the boundary of the j-th block the minimum and 
maximum limits of the effective transmissivities are defined 



rf Y := liminf rf Y := limsupr^. (10) 

71 >OQ — inn 



Hence, recalling that the Tj? Y, s are in nondecreasing order, 
we may notice that for any 5 > and for sufficiently large £ 



XY 



8<T$l 1)t+k <rf* +5, k = !,...£, (11) 



For each J, we are thus led to define two new sets of 
memoryless multi-mode lossy channels characterized, respec- 
tively, by the two sets of transmissivities {zf Y }j=i....J> and 



{Tj}j—i t ...j. Taking the limit £ — > oo while keeping J con- 



stant, their capacities can be computed as in Ref. 11311 yielding 



C 



7?> 



XY 



K :1 ) 



C 



3=1 



7?> 



XY 



Nj), 



3 = 1 



(12) 

where g(x) := (x + 1) log 2 (x + 1) — xlog 2 x I20I1 . The 
optimal photon numbers N_j and Nj are chosen in order to 
satisfy the energy constraint ( fTTT i and to guarantee the maxi- 
mum values of C_ and C respectively. Furthermore Eq. dTTT> 
shows that, one by one, each lossy channel entering the rhs of 
Eq. (0 can be lower or upper bounded by the corresponding 
channel of the two sets (this is a trivial consequence of the 
fact that a lossy channel can simulate those of smaller trans- 
missivity). Therefore the capacity of can be bounded by 
the capacities C and C of Eq. (fT2l . i.e. 



iff Y N 



3 J • 



(13) 



which applies for all J and for all XY. Taking the limit J — > 
00 and applying (0 we notice that the two bounds converge 
to the same quantity. Therefore we conclude that 



C = 



^.g(r(z)N(z)), 



(14) 



with N(z) being the optimal photon number distribution. Fol- 
lowing [13] it can be computed as N(z) = [t(z)(2 l/t ^ - 
where L is a Lagrange multiplier whose value is deter- 
mined by the implicit integral equation L n ^N(z) = N, 
which enforces the input energy constraint. In some limiting 
cases Eq. ( TBI admits a close analytical solution. For instance 
in the memoryless configuration e = 0, we get t(z) = rj, 
N(z) = N and thus correctly C = g(r]N) 11311 . Vice-versa 
for 77 = 1 (noiseless channel) or e = 1 (perfect memory chan- 
nel) we have t(z) — 1, N(z) = N and thus C = g(N) (per- 
fect transfer). Finally for r/ = (quantum shift channel) we 
get t(z) = e, N(z) = N and thus C = g(eN). For generic 
values of the parameters the resulting expression can be nu- 
merically evaluated, showing an increase of C for increasing 
memory e — see Fig. [2 a). 

Quantum capacity:- We proceed as in the previous case and 
use the results of Ref. lfl4ll for the quantum capacity on mem- 
oryless lossy channels to produce the following bounds on the 
quantum capacity of <&„ 

jt^ Y ^3)<Q<jf<l(Tf Y ^ (15) 



3=1 



J 



3=1 



max{0, g (tN) 
g((l — t)N)} is the maximal coherent information Il9ll 
the optimal photon number distributions Nj, Nj 



which holds for all J. Here q(r, N) 

and 

Nj can be com- 
puted as in Ref. Ill3ll . Finally we take the limit J — > 00 
applying Eq. (0 to the function q(r, N), yielding Q = 



4 




However generalizing it to non Gaussian inputs is problem- 
atic due to the infinite dimension of the associated Hilbert 
spaces 12 111 . Moreover, our results on the quantum capacity 



suggest that the channel is not forgetful if the input energy is 
unconstrained. 

We conclude by noticing that the optimal encoding strat- 
egy for the memory less channels which bound make use 
of coherent states 11311 . Since the latter are preserved by the 
encoding transformation Wa our results prove, as a byprod- 
uct, the optimality of coherent state encoding for the memory 
channel. 
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FIG. 2: (Color online.) (a) Contour plot of the classical capacity C as 
function of r\ and e for N = 8. (b) Contour plot of the unconstrained 
quantum capacity Qog, in the XY — EE, AE setups, as function 
of r\ and e (it diverges logarithmically at 77 = 1 and e = 1). (c) 
Distribution of the effective transmissivities r* Y for r\ — 0.7, e = 
0.3: the solid line shows the asymptotic distribution computed from 
Eq. (fSJ, the dots and the circles respectively show the distribution of 
transmissivities t^ e and t^ b for n = 20. 

Jo* If Q( T ( z )i N(z)), w i tn tn e optimal photon number dis- 
tribution N(z) to be computed numerically. 

A little thought leads to recognize that in the XY = 
EB, AB setups, where the output memory is assigned to Bob, 
there is at least one mode which is transmitted with unit effi- 
ciency for any value of e. In the case of unconstrained in- 
put energy this leads to infinite quantum capacity, implying 
that the limits n — > 00 and N — > 00 do not commute. A 
numerical evaluation indicates that, in the XY = EE, AE 
setups, the distribution of the transmissivities converges uni- 
formly to the function in (0. This implies that the formula (0 
can be applied even in the unconstrained case, yielding Qoo = 
Jo *" ^9(t(z)), where q(x) := max{0, log 2 x-log 2 (1 - a;)} 
- see Fig. EJb). 

Conclusions:- We have computed the capacities of a broad 
class of lossy bosonic memory channels without invoking 
their forgetfulness. Proving that the channel ([TJ is forgetful 
requires to show that in the limit n — > 00 the final state of 
the memory M (i.e. the state associated with the mode m' n ) is 
independent, in the sense specified in Ref. [2], on the memory 
initialization. A simple heuristic argument suggests that this 
is the case. The argument goes as follows: a photon entering 
from the input port mi of the setup has only an exponentially 
decay probability (erf) n of emerging from the m' n output port 
(this is the probability of passing through the sequence of n 
beam-splitters of of Fig. [T]>. Consequently the contribution of 
mi to the output state m' n is negligible for large values of n. 
If one restricts the analysis to Gaussian inputs with bounded 
energy this observation can be formalized in a rigorous proof. 



[1] C. H. Bennet and P. W. Shor, IEEE Trans. Inf. Th. 44, 2724 
(1998). 

[2] D. Kretschmann and R. F. Werner, Phys. Rev. A 72, 062323 
(2005). 

[3] K. Banaszek etal, Phys. Rev. Lett. 92, 257901 (2004). 

[4] E. Paladino etal, Phys. Rev. Lett. 88, 228304 (2002). 

[5] V. Giovannetti and S. Mancini, Phys. Rev. A 71, 062304 (2005). 

[6] M. B. Plenio and S. Virmani, Phys. Rev. Lett. 99, 120504 
(2007); New J. Phys. 10, 043032 (2008); D. Rossini, V. Gio- 
vannetti, and S. Montangero, ibid. 10, 115009 (2008). 

[7] V. Giovannetti and R. Fazio, Phys. Rev. A 71, 032314 (2005); 
D. Rossini, V. Giovannetti and R. Fazio, Int. J. of Quan- 
tum Info., 5 439 (2007); A. Bayat et al, Phys. Rev. A , 77, 
050306(R) (2008); V. Giovannetti, D. Burgarth, and S. Mancini, 
ibid. 79,012311 (2009). 

[8] V. Giovannetti, J. Phys. A 38, 10989 (2005). 

[9] N. Datta and T. C. Dorla s, J. Phys. A: Math. Theor. 40, 8147 
(2007); larXrv:0712.0722 !rquant-phl; A. D' Arrigo, G. Benenti, 
and G. Falci, New J. Phys. 9, 310 (2007). 
[10] G. Bowen, I. Devetak and S. Mancini, Phys. Rev. A 71, 034310 
(2005). 

[1 1] G. Bowen and S. Mancini, Phys. Rev. A 69, 012306 (2004). 
[12] Similar results can be also derived for their entangled assisted 
counterparts 

[13] V. Giovannetti et al., Phys. Rev. Lett. 92, 027902 (2004); V. 

Giovannetti et al, Phys. Rev. A 68, 062323 (2003). 
[14] M. M. Wolf, D. Perez-Garcfa and G. Giedke, Phys. Rev. Lett. 

98, 130501 (2007). 
[15] R. G. Gallager, Information Theory and Reliable Communica- 
tion (Wiley, 1968). 
[16] R. M. Gray, Toeplitz and Circulant Matrices: A Review, (Now 

Publishers, Norwell, Massachusetts, 2006). 
[17] A. S. Holevo, Probabilistic and Statistical Aspects of Quantum 

Theory (Amsterdam: North-Holland, 1982). 
[18] Va indicates the same transformation Vb when acting on the 

operators a k 's, i.e. V\a k VA = J2k' Okk' a h'- 
[19] A. S. Holevo and R. F. Werner, Phys. Rev. A 63, 032312 (2001). 
[20] The new sets can be treated as memoryless channels, as for n — > 

00 we have i = n/ J — > 00 copies of their transmissivities. 
[21] Also we notice that the notion of forgetfulness as introduced in 

|0] cannot be directly applied to the infinite dimensional case, 

and a proper extension would be needed. 



